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ON THE BEHAVIOR OF VELOCITY FLUCTUATIONS IN RAPIDLY ROTATING 

FLOWS 


S. S. GIRIMAJI* AND J. R. RISTORCELLP 

Abstract. The behavior of velocity fluctuations subjected to rapid rotation is examined. The rapid 
rotation considered is any arbitrary combination of two basic forms of rotation, reference frame rotation 
and mean flow rotation. It is recognized that the two types of rotating flows differ in the manner in which 
the fluctuating fields are advected. The first category is comprised of flows in rotating systems of which 
synoptic scale geophysical flows are a good example. In this class of flows the fluctuating velocity field 
advects and rotates with the mean flow. In the rapid rotation limit, the Taylor-Proudman theorem describes 
the behavior of this class of fluctuations. Velocity fluctuations that arc advected without rotation by the 
mean flow constitute the second category which includes vortical flows of aerodynamic interest. The Taylor- 
Proudman theorem is not pertinent to this class flows and a new result appropriate to this second category 
of fluctuations is derived, demonstrates that general fluctuating velocity fields are rendered two-dimensional 
and horizontally non- divergent in the limit of any large combination of reference frame rotation and mean- 
flow rotation. The concommitant ‘gcostrophic 5 balance of the momentum equation is, however, dependent 
upon the form of rapid rotation. It is also demonstrated that the evolution equations of a two-dimensional 
fluctuating velocity fields are frame-indifferent with any imposed mean-flow rotation. The analyses and 
results of this paper highlight many fundamental aspects of rotating flows and have important consequences 
for their turbulence closures in inertial and non-inertial frames. 

Keywords: Taylor-Proudman theorem, rotating turbulence, vortical flows 

Subject classification: Fluid Mechanics 

1. Introduction. Many important flows of geophysical and technological interest arc strongly influ- 
enced by rotation. It is generally believed that the fluctuating velocity fields in rapidly rotating flows 
are described by the Taylor-Proudman theorem (Greenspan, 1968, Tritton, 1977). The Taylor-Proudman 
theorem and one of its consequences, the geos trophic balance, have long been used successfully in many 
geophysical applications (Pcdlosky, 1982). 

Flow rotation can take many forms most basic of which are reference frame rotation and mean flow 
rotation (vorticity). In the former case, the flow occurs in a rotating reference frame, whereas, the latter 
category corresponds to an imposed mean flow rotation. A general rotating flow can be an arbitrary com- 
bination of these two basic flow types. In this article, the term rotating flow is used describe any flow that 
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is influenced by rotation; the terms reference frame rotation and mean flow rotation are the terms used to 
indicate the two basic types of rotation. The distinction between the two basic forms of rotating flow is best 
understood by considering one example of each flow: the homogeneous rotating shear flow and the elliptic 
streamline flow (Blaisdell and Sharif, 1996). These are two benchmark flows used to study the effect of 
rotation on turbulent fluctuations using direct numerical simulations (DNS). In the rotating homogeneous 
shear, the velocity fluctuations are subjected to a constant mean shear in a rotating reference frame. An 
elliptic flow is a homogeneous flow subjected to a constant mean strain and rotation in an inertial frame. 
The elliptic flow derives its name from the elliptic shape of the mean streamlines when the magnitude of 
mean rotation exceeds that of the mean strain. In the absence of strain, both the rotating and elliptic flows 
have circular streamlines. When the magnitude of strain and rotation are equal, the flow is a pure shear and 
the streamlines are rectilinear. 

Despite the apparent similarities between the flows with reference frame rotation and flows with mean 
flow rotation, the fluctuating velocity fields behave very differently, Blaisdell and Sharif (1996). The difference 
between the two flows increases as the rotation rate increases despite the fact that the mean streamlines in 
both cases become circular. An important kinematic distinction between the fluctuating fields in the two 
cases can be made in an inertial frame. In the rotating homogeneous shear the principal axes of the Reynolds 
stresses rotate at the rate of the mean rotation. In the elliptic flow, the principal axes of the Reynolds stresses 
are fixed. It is an explanation of this physical observation that motivates aspects of this article. 

The present state of development of the mechanics of rotating fluids does not distinguish between these 
two types of rotating flows. One important practical consequence of this fact is that the turbulence models 
calibrated in one class performs poorly in the other, Blaisdell and Sharif (1996). There is clearly a strong 
motivation for understanding the fundamental difference between the various classes of rotating flows and 
enhancing our ability to model the effect of rotation. 

The objective of this paper is to examine the behavior of various types of rotating flows in the limit 
of rapid rotation. The article is concerned with fundamentals; the mathematical distinction and physical 
interpretation of the various flow types. It is first recognized that the fundamental difference between the 
various flows lies, as is suggested by the physical observation alluded to above, in the manner in which the 
fluctuating fields are advected. Three types of fluctuations are identified. In the limit of rapid mean flow 
rotation, we establish for each type (i) the mathematical description and (ii) its distinctive kinematical and 
dynamical properties. The new results highlight some important fundamental aspects of rotating flows and 
have important mathematical and phenomenological implications for their turbulence modeling. 

In Section 2 the fluctuating energy equation is examined in the limit of rapid rotation to derive suitable 
approximations for the advection of fluctuating velocity and vorticity. In Section 3, the implications of these 
approximations on the fluctuation vorticity and momentum equation are derived for each flow class. The 
effect of reference-frame rotation and mean-flow rotation on two-dimensional fluctuations is examined in 
Section 4. Important turbulence modeling implications and physical insights are discussed in Section 5. We 
conclude with a summary of the results in Section 6. 

2. Rapid rotation: velocity fluctuation magnitude. We first demonstrate that the magnitude of 
the velocity fluctuations are, in the limit of rapid rotation, advected without change. This asymptotic result 
is then used to derive the permissible representations for the advection of the velocity and vorticity fields 
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in the rapid rotation limit. Establishing the form of the advection of the velocity immediately suggests a 
typology for rotating flows based on the rotation of the fluctuating velocity field or, as will be shown, on 
the rotation of the principal axes of the Reynolds stresses. This is a more unique characterization than the 
current distinctions of flow types based on the mean flow rotation or frame rotation. A classification of 
rotating flows into reference frame rotation and mean flow rotation is undesirable in as much as it involves 
the definition of the reference frame in the distinction. Moreover, such a classification cannot distinguish 
between the two flows when both arc considered in an inertial reference frame. 

The evolution equation for the fluctuating velocity implied by the Navier-Stokes equations in a non- 
inertial reference frame rotating at a rate is given by: 

(l) "(“ UjUtij + (Wij H - 2 f Tifcj + Sj, j XI j V 'Uiijj j 


where ^ represents the advection following a mean fluid particle: 

( 2 ) 


D 

Dt 


- + u^. 

dt 3 dxj 


Upper case letters represent mean quantities and lower case indicates fluctuations. The pressure in normalized 
by the constant density. The kinematic viscosity is given by v and e ijk represents the alternating tensor. 
Subscript ,j, represents partial derivative with respect to the j direction. Repeated indices imply summation. 
The mean strain rate, S*j, and the mean rotation rate, are given by 


( 3 ) 


S t j 


liUij + u^y, w i: = l -(u ld - 1 /,,,). 


The mean intrinsic vorticity Wij is defined as the sum of the system rotation and mean flow rotation: 


(4) Wjj — -f 2 

The equation for the magnitude of the velocity fluctuations is obtained from equation (1): 

( 5 ) + 1 = ~Wij u j u i ~ SijUjUi - p,iUi + VUiUi,jj . 

= -SijUjUt -p,iUi + I 'UiUi,jj . 


This equation will be used to show that in the limit of rapid rotation that the magnitude of the velocity 
fluctuations are approximately constant. 


2.1. The rapid rotation limit. We are interested in the evolution rate of the magnitude of the velocity 
fluctuations in the limit of rapid rotation or, more precisely, large intrinsic vorticity. Throughout the paper, 
as in the case of the Taylor-Proudman theorem, the rapid rotation limit implies: 


( 6 ) 


W 


« i; 


w <<1 ’ 


where, 5, W and oj represent, respectively, the norms of the mean strain rate, intrinsic mean vorticity 
and fluctuating vorticity. The square of the fluctuating vorticity, u> 2 can also be interpreted as the energy 
cascade rate. The large intrinsic vorticity analysis includes, of course, that rapid reference frame rotation 
as well as the large mean flow rotation cases. Our analyses are to be understood, in the same sense as the 
Taylor-Proudman theorem, as an asymptotic analysis in the limit specified. 
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On normalizing by the intrinsic mean vorticity, the equation for the magnitude of the velocity fluctuations 
becomes 


(7) 2W^ U j( U i U i)ij ] ~ yy[Sij u j u i (P u i),i + V u i u i)jj ]• 

The terms on the right hand side of the above correspond to the production, pressure transport, and viscous 
dissipation of the energy of the fluctuations. The production and dissipation will scale with the cascade rate, 
which in the rapid rotation limit, is small (equation 6). The turbulent transport terms do not contribute 
towards the overall growth of fluctuating energy as they merely transport energy from locations of production 
to areas of dissipation. If production and dissipation are individually small, these effects must be small. Thus, 
on the time scale of the intrinsic vorticity 


( 8 ) 


1 

W 


UjUiJ 


0; yy SlJ M, Uj 


0; ^v VUi Ul ' n 


0; w 


0. 


Thus, on the time scale of the intrinsic vorticity the magnitude of fluctuating velocity following a fluid 
particle changes slowly: 


( 9 ) 


1 D f ^ 

wWt {n ' Ui) 


0. 


That is, the magnitude of the velocity is approximately constant. Given that the magnitude of the velocity 
is approximately constant, the advection must be of the general form 


( 10 ) 



'R% j U'j i 


where 7 Zij is a general antisymmetric tensor so that, in forming the energy equation, TZijUjUi = 0. The form 
of (10) is, of course, recognized as the pure rotation of a vector field. 

The arbitrary rotation tensor 7 Zi 3 can be decomposed into mean and fluctuating : 


(ii) 


R'ij — Rij ~b Fiji Aj 


where Ri 3 represents a systematic (mean) rotation of the velocity field while represents random rotations 
about the mean rotation by the fluctuating field. The mean rotation rate can be expected to scale with the 
mean intrinsic vorticity and the random rotation with the fluctuating vorticity: 


( 12 ) 


Jr|_ 

1*1 


H«i 

W ’ 


as is consistent with the inequalities already assumed. Therefore, on the timescale of intrinsic vorticity, 


(13) — Ui « RijUj. 

This equation, describing the pure rotation of the velocity field at some as of yet unspecified rotation rate, 
is the fundamental equation from which we shall obtain several results pertinent to flows undergoing rapid 
rotation. The rotation rate Ri 3 will be called the kinematic rotation tensor. 

The implication of equation (13) for the rotation of Reynolds stresses is straightforward: 


(14) 


— Rik{UjUk) H - Rjk • 
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The kinematic rotation rate of the fluctuating velocity field is, on the time scale of the rotation, equal to 
the rotation rate of the principal axes of the Reynolds stresses. Thus, statistically R XJ can be understood 
or even called the rotation rate of the Reynolds stresses. 

The implication of equation (13) for the fluctuating vorticity is easily obtained by taking its curl: 

(15) ~Dt^ i ~ 2 ^ J ~ Wijujj = — RjUij). 

Note that the advection of the vorticity is not a simple rotation: there are very complicated production terms 
involving the fluctuating velocity gradients. In the above equation, the following relationships between an 
antisymmetric tensor and its associated axial vector have been used: 

(16) Wi= e, jk W kj - Wij= - ±e ijk W k ; 

1 

OJi — EijfaLdkj J Uij — 2 

Ri = EijfcRhj] Rij — — EijkRk- 

To summarize: In the limit of rapid rotation we have shown that the advection of the velocity approaches 
a pure rotation. There are two important implications of this result. The first is that the rotation of the 
principal axes of the Reynolds stresses is a very close approximation to the rotation rate of the velocity field. 
Secondly an expression for the advection of the vorticity, consistent with a pure rotation of the velocity 
field, can be obtained. The advection of vorticity is not a simple rotation: there are source terms that may 
contribute to the growth of enstrophy. Equations (13) and (15) form the foundation of this paper, as they 
will be used in the next section in the momentum and vorticity equations to infer important properties of 
the fluctuating fields in the limit of rapid rotation. 


2.2. Types of kinematic rotations. In the limit of rapid rotation it was shown that the advection 
of the velocity field is given by a pure rotation, 

D 


(17) 


Dt 


RijUj 


where Rij, by definition, is the rotation rate of the Reynolds stress tensor. We will now consider the values 
of Rij associated with reference frame rotation, mean- flow rotation and arbitrary rotation. The value of 
Rij will depend upon the reference frame of observation. In the inertial frame, it will be seen that Rij = 0 
corresponds to mean-flow rotation and R^j = W X j corresponds to reference frame rotation. Knowing the 
kinematic rotation rate tensor, we obtain expressions for the advective derivatives of velocity and vorticity 
for the various flow classes. These expressions will be used to approximate the advective derivatives in the 
rapid rotation analysis of the vorticity and momentum equations. 

Reference-frame rotation . It is known from observations of rotating shear flows that the Reynolds stress 
axes of this type of flow are fixed in a rotating reference frame. Therefore, 


(18) 


Rij — 0 


{ Dt 


Ui]R 


°. [§- t Ui]R 


0. 


In an inertial frame, the Reynolds stress tensor rotates with the mean flow at the rate of mean vorticity 
If Wij is the mean flow rotation tensor in an inertial frame then 


(19) 


Rij = 


D 


D 


£)j. Ui ~ WijUj, ~ WijWj — 2 WjUJjj. 
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As these fluctuations are Advected With the Mean flow Rotation they will be denoted AWMR fluctuations. 

The mean Oldroyd derivative is a form of the substantial derivative convecting and rotating with the 
mean flow. The advective derivative in a rotating reference frame is exactly equal to the Oldroyd derivative 
in an inertial frame. For the AWMR fluctuations, the advective derivative is negligible in the rotating frame 
and the Oldroyd derivative is negligible in the inertial frame. 

Mean- flow rotation . In the elliptic mean flows already mentioned, it has been observed that the principal 
directions of the Reynolds stress are fixed in an inertial frame. The rotation rate of the Reynolds stress is 
zero, thus 

( 20 ) Rij =0 => ~ ~Dt^ i ^ = 2 WjUjj. 

As these fluctuations are Advected With Out the Mean flow Rotation they will be denoted AWOR fluctua- 
tions. 

General mean flows. Fluctuations in the presence of arbitrary combinations of reference frame rotation 
and mean flow rotation are of substantial interest. The rate of rotation of the Reynolds stresses is not, 
a priori known. Let the arbitrary rotation rate of the Reynolds stresses in a given frame of reference be 
denoted by Rfj. In that frame of reference, the fluctuating field description is 


( 21 ) 


Rij — R t j ; 


D 


&ij u V 


Dt 


Rj)uij + WijOjj = -(WjUj t i - RjUij). 


This class of fluctuations will be called AWAR fluctuations since it represents Advection With Arbitrary 
Rotation. 

Examples of different classes. The specification of a mean flow rotation does not constitute sufficient 
information to uniquely determine, a priori, the type of fluctuating velocity field. However, a classification 
of flow type is possible in many important flows based on experimental/numerical observation. Synoptic 
scale geophysical flows in the atmosphere are very visible examples of AWMR fluctuations. All flows in 
rotating systems, including many technological processes crystal growth, centrifugal separation, rotating 
turbomachinery — are likely to fall into the AWMR category. In contrast, flows in stationary systems subject 
to mean flow rotation may belong to the AWOR category. Examples of such flows are many of the vortical 
aerodynamic flows — the trailing vortex, the flap-edge vortex, leading-edge vortex — and elliptic flows 
(Blaisdell and Shariff, 1996). Examples of AWAR category include elliptic flows in rotating reference frames, 
cylindrical Couette flow with the two cylinders rotating at different speeds and a cyclone which rotates with 
its characteristic angular speed on the rotating frame of the earth. In these flows, both reference- frame and 
mean- flow rotation are important. 


3. Rapid rotation limit: vorticity and momentum equations. In the previous section a possible 
typology for rotating flows according to the value of the rotation tensor, Rij, in the ansatz for the advective 
derivative of the velocity field, was indicated. This naturally suggested three different types of fluctuations 
which we have labelled AWMR, AWOR and AWAR. In this section the properties of AWMR, AWOR and 
AWAR fluctuations are, without regard to their physical origin, investigated. The kinematic ansatz for the 
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advective derivative will now be used, for the first time, in the dynamical equations. It should be emphasized 
that the he objective of this section is not an a priori categorization of a particular rotating flow, but rather 
an a posteriori examination of important properties of the diverse fluctuating fields. The implications for 
the dynamical equations of diverse Rij are explored. 

3.1. AWMR fluctuations. In this subsection AWMR fluctuations are investigated. The analysis for 
the AWMR fluctuations is, in fact, a derivation of the well known Taylor- Proudman theorem, Proudman 
(1916), Taylor (1917), Greenspan (1968), TYitton (1977). This will serve as a familiar introduction and 
set precedent for procedures that will be used to obtain results in subsequent developments. The Taylor- 
Proudman result will be derived in the familiar rotating frame here and in an inertial frame in Section 
3.3. 

The incompressible Navier- Stokes equations in a non-inertial reference frame rotating with an angular 
velocity $1* is given by 


( 22 ) 


tq + u j^t fj + ^ e ikp^k U p — P U VUi yjj . 


The instantaneous values of velocity and pressure are represented by u* and p* respectively. In this section 
on/y, asteriks will be used to represent dimensional quantities. Dimensionless variables are represented by 
the symbols without asteriks. 

The equations are made nondimensional with length and velocity scales i and u c respectively. The time 
scale r c is formed from i and u c . The magnitude of the reference frame rotation rate is given by fZ c . The 
Navier-Stokes equations become: 


(23) 


[ ^ [R°] 'U'jUiij + 2 £ikp£lkUp — p,i + [Ek]Uiyjj . 

TAlr WriZr 


Here the Ekman number, Ek — and the Rossby number, Ro = represent, respectively, the relative 
importance of viscous and inertial forces to Coriolis forces. It is also useful to view the Rossby number as 
an indication of the magnitude of flow vorticity relative to the background vorticity due to frame rotation. 

Rapid rotation limit of the vorticity equations .. Taking the curl of the Navier-Stokes equations one 
obtains the vorticity equation 


(24) 


— —] + [7?o] Ujtx>i,j = \Ro] UjUi,j -I- 2f} J u I , J + [Ek] w*,*,- , 

T c il c ’ 


where ujk = £kiqUi, q . In the limit of rapid rotation (ft c — » oo) the Rossby and Ekman numbers vanish and 
the fluctuating vorticity equation implied by the Navier-Stokes equation becomes 


(25) 




u>i , j . 


The advection of the fluctuating vorticity for this class of fluctuations can be approximated to zero as given 
in equation (18). This approximation is valid for vorticity fluctuations whose timescale r c is much larger 
than the timescale of the rotation of reference frame (1/£"2 C ), leading to r c £l c — ► oo and 


(26) 


2 CljUi,j — ► 0 . 


This is the Taylor-Proudman theorem and it is a statement valid only for those vorticity and velocity fluctu- 
ations that evolve slowly with respect to the frame rotation rate. It is valuable to consider the geophysical 
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point of view. A flow structure consistent with 2ttjUi,j — * 0 is necessary to keep, as is observed in nature, 
the relative flow vorticity (Rossby number) small. 

The Taylor-Proudman result has several important consequences (Lighthill, 1966, Greenspan, 1968, 
TYitton, 1977). Without loss of generality assume that the axis of the frame rotation coincides with the ‘3 1 
coordinate axis. The Taylor-Proudman result indicates that 


(27) O jUi,j —> 0 => Ui , 3 = 0 : 

the fluctuating field is invariant along the axis of rotation. The fluctuating velocity field is two-dimensional 
Ui = Ui(xi,x 2 ). Since 1 * 3,3 = 0, the continuity equation yields 


(28) 


Uui + 1 * 2,2 = 0 . 


The flow is nondivergent in planes perpendicular to Slj, 

Rapid rotation limit of the momentum equations. In the rapid frame rotation limit, the momentum 
equation reduces to 


(29) 


_ ]ui + 2 EikpQk^p ~~ Pii ■ 

• j l p U c j Z £ 


The advection of the velocity can be approximated as zero; for AWMR fluctuations Rij = 0 in the rotating 
frame equation (18). More precisely, for fluctuations that evolve slowly in the rotating frame, r c Q c — > 00 
and [ 777 ] iii — * 0 , leading to 

(30) 2 EikptokUp = p,i - 

"C 


The scales of the fluctuating motion that evolve slowly in the rotating frame satisfy, to leading order, the 
geostrophic balance — a balance between Coriolis and pressure forces. Using the Poisson equation for 
pressure the pressure can be shown to scale as p* ~ u c Q c £. As a consequence of the geostrophic balance the 
flow is described by a streamfunction. Multiplying the geostrophic balance by £i q ifli and contracting leads, 
in dimensional variables, to 


(31) 


1 n; . . 

U P- 2 £iqk p ,i+U k . 


The fluctuating AMWR velocity field is, to leading order, described by a streamfunction: 


(32) Up — £pqk&k' l fiq “b 'iptp — ^qpk^k^lq "b Qp^qlpiq 

where ip = ip (x, y). The streamfunction and pressure field arc proportional; isobars and streamlines are 
parallel. Additional details can be found in Ristorcelli (1997). 

In the present rotating coordinate system, AWOR fluctuations undergo a kinematic rotation. Therefore, 
AWOR fluctuations evolve on a time scale of the order of leading to r c n c — ► 1 . This is not consistent 
with r c fi c — ► 00 as required for the Taylor-Proudman theorem. 
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3.2. AWOR fluctuations. The AWOR fluctuations are now investigated. This is the type of fluctu- 
ations in which the velocity is advected but not rotated by the mean flow; Rij = 0. These fluctuations are 
most conveniently analyzed in an inertial reference frame. The momentum equation for the fluctuations is 
given by 


(33) 


Du* 

Dt 


■<+Kis: k +w* k ]= - P ‘, t +uu* 


r 33 > 


where indicates derivative following a mean fluid particle. The similarity with the momentum equations 
in the rotating frame (22) and (33) is clear: the mean vorticity now plays the role of the Coriolis force. It is 
important to note that the mean vorticity is twice the angular rotation: W£ = 2 fl*. 

Rapid rotation limit of the vorticity equations . Taking the curl of the momentum equation produces the 
equation for fluctuating vorticity: 


(34) 


Du>* „ * 

-dT + u W' 


x i ->3 


+ W J u i ij + (&ij + W*j) + v u* ,jj , 


Normalizing the fluctuating vorticity equation with characteristic length, time, velocity scales, and charac- 
teristic mean rotation rate r c , u c , W c produces 


(35) 


l r yy ] Dt “ WjUilj -h WijUJj + [ w ] SijLJj + 


[Ro w ]ujjUi,j + [Eky,] C0i 


133 * 


The corresponding Ekman and Rossby numbers are Ek — i and R° = wJ * Here ~ Wij R'o and 
= SijSij. The ratio [S c /W c ] determines the ellipticity of the flow. In the limit rapid mean rotation or 
large mean vorticity, W c — > oo, the Ekman, Rossby and ellipticity parameters vanish and 


(36) 


, Du>i 


IT) Ui "hJVijUJj 


t c W c Dt 

If the AWOR ansatz , equation (20), for the advective derivative of the vorticity fluctuations is applied, 
\WjUj y i one obtains the result 


p 

Dt 


(37) 


WjUiij 0. 


This is equivalent to the Taylor-Proudman result with Wj replacing fij. The velocity field is again two- 
dimensional being independent of the coordinate along the axis of the mean vorticity. The flow is also 
nondivergent in planes perpendicular to the rotation. 

Rapid rotation limit of the momentum equations . While both AWMR and AWOR fluctuations are both 
independent of the coordinate along the axis of rotation and are both horizontally nondivergent, they arc 
very different in other important aspects. This difference is seen in the form of the advective derivative 
which manifests itself dynamically in a different “geostrophic” balance of the momentum equations. 

The momentum equations for the AWOR fluctuations, in the large mean vorticity limit, arc 

(38) [^rl §^ + u kWik = - Pti . 

If the AWOR ansatz , equation (20), for the advective derivative of the velocity fluctuations is made one 
obtains the following balance 

(39) u k W ik = ^e ikq W q u k w - p,i - 


9 



Contracting the above balance with £ij q W q yields the streamfunction for the AWOR fluctuations: 


(40) 


XLj — 2, £{j q W q “b U q \fy • 


which in dimensional units becomes 


(41) 


W* 

2e ' kq w? p 


+ < 


w? ’ 


Pressure plays the role of a streamfunction as it did in the geophysical case. However, there is a factor of 
four difference in the prefactor of pressure in the two cases, equations (32) and (41). This reflects important 
dynamical differences between AWMR and AWOR fluctuations. 

Consider the AWMR fluctuations of §3.1 in an inertial frame. Due to their rotation, they change rapidly 
on a timescale, r c , which is small r c ~ For AWMR fluctuations, r c W c — + 1 and the analysis of the 
present section is not possible as r c W c — ► oc is required. Thus AWMR and AWOR fluctuations, despite 
the same implications for the vorticity equation and a “twodimensionalization” of the flow are distinct and 
mutually exclusive classes of velocity fields. 


3.3. AWAR fluctuations. The behavior of AWAR fluctuations in an inertial reference frame is now 
considered. The fluctuating velocity field evolves according to 


(42) 


~-u * + u*u* 




r 4 - W*U* + <?*■ 


' P n H - v U± ijj 


In the limit of rapid rotation, the effects of mean strain rate, viscosity and nonlinearity can be neglected 
leading to 


(43) 


D 

rn 


W*u* 


The fluctuating vorticity equation in the rapid mean rotation limit is 

D 


(44) 


Dt 


Ui 




These equations are now examined by invoking the approximations for the advective derivatives given in 
equation (21). 

Rapid rotation limit of the vorticity equations . The advection of the AWAR vorticity fluctuations is 
approximated as, equation (21): 


(45) 




When this simplification is substituted into equation (44) we obtain: 


(46) + o, 

for all combinations of W* and If W* and are coincident with the 3-axis, then the fluctuating velocity 
field is two-dimensional and horizontally non-divergent: 


( 47 ) = 0 ; — 0 , 

as is the case for the Taylor-Proudman relevant to AWMR fluctuations. 
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Rapid rotation limit of the momentum equations. The advective derivative of the AWAR fluctuating 
velocity can be approximated as, equation ( 21 ), 


(48) 

From equation (43) one obtains 


Tn< " 


(49) 


(J?“ + 


This is a “geostrop hie” balance for AWAR fluctuations. Such balance of the momentum equations, as shown 
in §3.1 and §3.2, leads directly to a streamfunction representation for the components of the velocity in 
planes perpendicular to the axis of rotation. The streamfunction representation is, of course, due to the 
horizontal nondivergence of the flow. 

The results for AWOR fluctuation can be obtained from above results by recognizing that for AWOR 
fluctuations Rfj = 0. The results for the AWMR fluctuations, the Taylor-Proudman theorem in an inertial, 
is now derived. 

AWMR fluctuations in an inertial frame. The advection of these fluctuations in the inertial frame is 
given by (equation 19): 


(50) 


— u*= W*u*- 
Dt ' ** 


§i“' - W M- 


Here Wf- is the mean- flow rotation rate. This indicates that the Oldroyd derivative of these fluctuations are 
small in the inertial frame, consistent with the advective Lagrangian derivative being small in the rotating 
frame. By recognizing R = W*j in equations (46) and (49), we obtain the Taylor-Proudman theorem and 
its corresponding geostrophic balance in an inertial reference frame: 


(51) 




0 ; 




e ipj w; u * 


- p,i ■ 


That this result is identical to the Taylor-Proudman geostrophic balance equation (30) is seen readily by 
recognizing that W£ — 2 

4. Discussion. Some useful insights into the results derived arc now given. Further, the implications 
of the results for developing turbulence models for rotating flows is also discussed. 

Kinematics of the fluctuating vorticity field. The fluctuating vorticity equation indicates that all classes 
of fluctuating velocities are rendered two-dimensional in the rapid rotation limit (equations 47, 46). The 
two- dimensionality of the fluctuating velocity field has important implications. If the mean vorticity vector 
is oriented along the 3- axis, then the fluctuating velocity field is such that u ^ 3 = 0. Thus W$Ui$ ~ 0 for 
all fluctuations in the limit of rapid rotation. When this result is introduced into the vorticity evolution 
equation (44), we obtain 

(52) -^w t = + UjWij ss uijWij. 

The fluctuating vorticity field of a two-dimensional flow undergoes a kinematic rotation in the limit of rapid 
rotation. The rotation rate of vorticity is equal to the mean flow rotation rate. In associated enstrophy 
equation is 

(53) = OJiUfjWij = 0. 
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The last equality follows from the antisymmetry of Wij. On the timescale rotation, the enstrophy is constant. 
If the fluctuating vorticity is small in comparison to the mean intrinsic vorticity at the initial time, it will 
continue to be small on the timescale of the intrinsic vorticity. Thus, a fluctuating velocity field that satisfies 
the rapid rotation limit assumptions specified in equation (6) will preserve the validity of these assumptions 
for all time. 

The pressure-rotation balance. In a general rotating flow, the geostrophic balance between the pressure 
and rotation- related forces depends upon the kinematic rotation rate of the fluctuating velocity field (equation 
49). Consider AWMR and AWOR fluctuations in an inertial frame, subject to identical mean- flow vorticity 


W*j. The fluctuating pressure gradients for the two cases are: 


(54) 


AWMR: 

AWOR: 


P,i ~ T\V{j Uj 5 
P,i ^ j tlj . 


For a given intrinsic mean-vorticity, it takes twice as large a (fluctuating) pressure force to sustain AWMR 
velocity fluctuations, as it does to maintain AWOR velocity fluctuations. The additional force is required in 
the AWMR case to rotate the orientation of the velocity field. 

Two-dimensional fluctuations . One well known result in the mechanics of rotating fluids, (Taylor, 1917, 
Lighthill, 1966, Hide, 1977) is that two-dimensional velocity fields are frame- indifferent. Any two-dimensional 
fluctuating velocity field, U{ = U{(x 1 , 22 ), ^at is a solution of the Navier-Stokes equations in an inertial 
system is also a solution of the Navier-Stokes equations in a rotating frame for any arbitrary (constant) 
rotation rate. The pressure fields for the two evolutions are not the same but proportional by a constant 
scale factor. 

In the previous section it has been seen that in the limit of rapid rotation AWMR, AWOR and AWAR 
fluctuations become independent of the coordinate along the axis of rotation - two-dimensional. Thus 
AWMR, AWOR and AWAR fields are frame indifferent. The shortest proof of the material frame indifference 
of two-dimensional fields uses the vorticity equation. The exact fluctuating vorticity equation for arbitrary 
(small or large) frame or mean flow rotation rate, modulo viscosity, is 


(55) 


D 

— UjU>i,j= [Uj + Wj\u i , j -ft OjWij. 


If the flow is two-dimensional independent of coordinate along the axis of rotation then = (Wj + 

2fij)uj,j = 0 and the vorticity equation becomes 

D 


(56) 


Dt 


U>i + Uj(jJ{,j — U)jUi,j *f WjjUJj. 


The vorticity equation and the kinematic Biot-Savart relation V 2 u* = £ijkWk,j are all that are required to 
evolve the flow. The vorticity equation and Biot-Savart are both independent of the rotation the reference 
frame, Qj. The flow is frame-indifferent in the presence of mean vorticity. 

A two-dimensional fluctuating velocity field is more easily evolved by solving the vertical vorticity equa- 
tion. Consider the evolution of the 3-component of vorticity. Noting that W 3 j = 0 and that UjU 3 ,j ^ 0 for 
velocity fields described by a streamfunetion one obtains 

D 


(57) 


Dt 


LJ3 -f UjU 3 ,j 


0 . 
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As (Ji, 3= 0 the equation is independent of the vertical velocity, U3. Thus 3 = 0 is a closed nonlinear 
equation for the streamfunction from which the horizontal components, U\>U 2 , of the velocity are fully 
determined. Recall that U 3 = — V 2 rp as is readily obtained from the curl of the velocity once the velocity is 
expressed in terms of the streamfunction. 

The third component of the velocity field, U3, is obtained from its momentum equation. With ui l3 = 0 
and p, 3 = 0 the 3 momentum equation devolves into the linear equation describing the advection of a passive 
scalar. The passive scalar is u 3 . For arbitrary frame rotation, u 3 satisfies 

(58) ~^u 3 + ui« 3j i + 1 ^ 3 , 2 = 0. 

This u 3 equation is, of course, independent of the mean frame rotation. 

The fact that two-dimensional velocity fields are materially frame indifferent, in the context of rotat- 
ing flow systems, has been known for some time (Taylor, 1917, Hide, 1977). The utility of the principle 
of materially frame indifference in the limit of two-dimensional flows has been suggested as a type of re- 
alizability principle for turbulence models by Speziale (1981). Additional amplifications on the topic of 
two-dimensional frame indifference are given in Speziale (1989, 1997), Kassinos and Reynolds (1994), and 
Ristorcelli (1997). The frame invariance of the modeled second moment equations has been implemented for 
single-point turbulence closures relevant to AWMR flows of the geophysical type in Ristorcelli, Lumley and 
Abid (1995). 

Implications for Reynolds stress closures. In the Reynolds stress closure modeling of turbulent flows, one 
requires a closure representation for the pressure- strain covariance, Launder, Reece and Rodi (1975). The 
present results have important implications for pressure-strain closures. These results describe the behavior 
of the Navier-Stokes equations in the limit of rapid rotation. The Reynolds stress equations, which are 
derived from the Navier Stokes equations, must, in the limit of rapid rotation, be consistent with these 
results. 

In Ristorcelli et al. (1995) the geostrophic balance, a diagnostic relationship between pressure and 
velocity for AWMR flows for which Taylor-Proudman is relevant, 

(59) 2 WijUj » -p,*; 

has been employed to develop a sophisticated representation for the pressure-strain. This model has demon- 
strated advantages over models not consistent with Taylor-Proudman results in rotating flow calculations. 

One of the major results of this article is a new rotation-pressure balance of the momentum equations 
for general rotating flows: 

(60) (Wij ~b R-)Uj ~ ~ Pn ‘ 

This result is used by Girimaji (1997) to develop a pressure- strain covariance model for mean vorticity 
dominated flows. The resulting pressure- strain model is the first of its class to distinguish between rotating 
shear and elliptic flows, and produces excellent qualitative and quantitative agreement with DNS data of 
elliptic and rotating flows. 

Estimation of Rf 3 . Many important kinematic and dynamic properties of the fluctuating fields in a 
rapidly rotated mean flow can be obtained provided that the Reynolds stress rotation rate (Rfj), which is 
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also the kinematic rotation rate of the fluctuating field is known. It would be desirable to specify in terms 
of mean velocity field variables. The kinematic rotation rate is likely to depend upon many factors: initial 
condition, boundary condition and production- mechanism of the fluctuations. For a homogeneous rotating 
flow in a fully developed state, the kinematic rotation rate is likely to depend only on the fluctuation 
production- mechanism, which is the mean strain rate. We now investigate the dependence of Rfj on the 
mean strain rate. This qualitative analysis will provide a better physical insight into the difference between 
the fluctuating fields of rotating shear flows and elliptic streamline flows. 

From observations, it is well known that the angle between the principal axes of the mean strain and the 
Reynolds stresses is constant in equilibrium turbulence. The Reynolds stress anisotropy attaining a constant 
equilibrium value is a manifestation of the unchanging orientation between strain rate and Reynolds stresses. 
The equilibrium behavior is observed in flows with reference frame rotation, as well as those with mean flow 
vorticity. The rate of rotation of the Reynolds stresses must, then, coincide with the rotation of the strain 
rate: 


(61) 




Er 

Dt %3 ' 


• In elliptic streamline flows, the strain rate is constant in magnitude and orientation. Thus in an 
inertial reference frame: 


(62) R% = 0. 

This indicates that velocity fluctuations in the elliptic mean streamline flow belong to the AWOR 
class. 

• In rotating shear flows, the strain rate eigenvectors rotate with the mean flow when observed in an 
inertial reference frame. The intrinsic mean vorticity, Wij , and the rotation rate of the strain rate 
are equal. 

(63) R% = W tj = 2 EikjSlj 

where Qj is the reference frame rotation rate. Therefore, the velocity fluctuations in the rotating 
shear are likely to belong to the AWMR category. 

5. Summary and Conclusion. The behavior of fluctuating motions subject to a general rapid rotation 
has been investigated. The forms of rotation considered include reference- frame rotation, mean-flow rotation 
and any arbitrary combinations of the two basic forms. 

In the rapid rotation limit it has been demonstrated, §2.1, that the fluctuating velocity field advection, 
to leading order, can be approximated by a kinematic rotation given by 

D 

(64) Dt Ui ~ RijUj ‘ 

The kinematic rotation rate tensor Ri 3 was shown to be equal to the rotation rate of the Reynolds stress 
tensor. The advection of fluctuating vorticity corresponding to a kinematically rotated velocity field is given 

by 

(65) ~ 2^^^ _ WijUtj = — RjUij). 
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These two equations form the foundation of this study. The introduction of these kinematic approximations 
into momentum and vorticity equations indicated the behavior of the fluctuating velocity and pressure fields 
in the rapid rotation limit. The rotating flows were characterized according to the rotation of the 
Reynolds stress tensor, and the behavior of each class was investigated in §3. 

The results of the analysis for each class are summarized below: 

• Rij = Wiji This fluctuating velocity field is advected with mean-flow rotation rate (AWMR). The 
AWMR class of fluctuations can be associated with reference-frame rotation flows. The vorticity 
and momentum equations, in this case, lead to the classical Taylor-Proudman theorem and the 
geostrophic balance, respectively. The Taylor-Proudman theorem results from assuming that the 
Lagrangian derivative, in the rotating reference frame, is small: in the rotating frame, on the time 
scale of the rotation, the flow is frozen. This is tantamount to assuming, in an inertial frame, that 
the Oldroyd derivative is small. In an inertial system, these fluctuations while rapidly varying, 
are merely kinematically rotated. In the limit of rapid rotation the fluctuating velocity field is 
two-dimensional independent of the axis of rotation - and horizontally non- divergent. The Taylor- 
Proudman theorem and its consequences have substantial empirical validation: these results arc 
essential to understanding large (synoptic) scale geophysical flows taking place in the atmosphere 
and the ocean. The geostrophic balance is used to produce weather maps. 

• — 0: This fluctuating velocity field is advected without rotation by the mean flow (AWOR). 
This class of fluctuations are likely to be associated with flows which have mean- flow rotation and 
no reference frame rotation. For AWOR fluctuating velocity field the Lagrangian derivative, in an 
inertial frame, is small; in the inertial frame, on the time scale of the rotation, the flow is frozen. 
In the limit of rapid rotation, the vorticity equation indicates that the fluctuating velocity field is, 
again, two-dimensional and horizontally non- divergent. However, the momentum equation indicates 
a different balance between rotation and pressure forces. These results appear relevant to turbulence 
closures of highly vortical flows of engineering interest. 

• Rij = R^j : This class of fluctuations will be encountered in flows which are subjected to both 
reference- frame and mean-flow rotations. The fluctuating velocity field is, again, two-dimensional 
and horizontally non- divergent. The pressure-rotation balance, however, depends on the proportion 
of reference- frame and mean-flow rotation. 

The asymptotic limits of rapid rotation obtained in this paper are of substantial fundamental interest as 
can be seen from their physical meaning, as well as the kinematic and dynamic consequences. These asymp- 
totic states of rapidly rotating flows have important implications for turbulence modeling Speziale (1981, 
1989, 1997), Kassinos and Reynolds (1994), Ristorcelli et al (1995), Ristorcelli (1997). To be consistent 
with Navier-Stokes equations, pressure-strain models of rapidly rotating flows must satisfy the geostrophic 
constraint as implied by the Taylor-Proudman theorem. Ristorcelli et al (1995) have employed this con- 
straint to develop a pressure- strain covariance model for flows in rotating systems. In contradistinction, in 
flows subjected to large mean vorticity, pressure-strains model must be consistent with a new geostrophic 
constraint derived in this paper. This constraint is used by Girimaji (1997) to derive a pressure- strain model 
able to predict the peculiar behavior of Reynolds stresses observed in the important elliptic streamline flows. 
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